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^ \ Possible non-standard top-quark interactions with the gluon and the W boson 

induced by 577(3) x SU(2) x J7(l) gauge-invariant dimension-6 effective operators 
are studied for hadron-collider experiments. Current limits on top-gluon couplings 
are presented by using the latest experimental data of tt productions at the Teva- 
tron and the Large Hadron Collider (LHC). The optimal-observable procedure is 
applied to the charged-lepton distributions in pp — > ttX — » £ + X' {t = e or u.) at the 
LHC in order to estimate the expected statistical uncertainties in measurements of 
those non-standard top-gluon and top-VF couplings that contribute to this process 
in the leading order. 
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1. Introduction 

The Large Hadron Collider (LHC) is currently one of the most powerful facilities 
for probing possible new physics beyond the standard model (BSM) pp. In collider 
experiments such as the LHC, there could be two kinds of approaches to new- 
physics studies: One is a direct search for new particles which do not belong to the 
standard model (SM). The other is an indirect search, which tries to find certain 
deviations from the SM predictions in processes having only well-known particles 
in their initial and final states. 

In the former case, any non-standard particles once being discovered, we would 
get a lot of information about new models which contain those particles. However, 
it is not that easy to enable such high-energy collisions which are enough energetic 
to create new heavy particles. On the other hand, in the latter case, observable 
deviations might appear even at an energy far below a typical new-physics scale 
since they could be induced by quantum effects of non-SM heavy particles. Of 
course a large number of events are needed in order to establish a new-physics 
evidence at a statistically-significant level, because such signals would generally be 
quite small, considering the remarkable success of the SM in various phenomena 
below O(10 2 )-GeV scale. 

Standard particles have been satisfactorily "re" discovered today except for the 
Higgs boson at the LHC, including top quarks. This is the first time for the 
top quark to be measured at a facility other than the Tevatron although almost 
twenty years have passed since its discovery there [2j [3] . This situation, in which 
more precise tests of top-quark properties can be continued, is quite meaningful for 
new-physics studies after the Tevatron is shutdown: The top quark is the heaviest 
particle of all standard particles, indicating that precise measurements of this quark 
might be able to open windows for the BSM j3]. 

In various direct searches of new particles up to today, any plausible candidates 
have not been discovered yet, which definitely requires indirect BSM searches. Con- 
sidering this circumstance, we pay attention in this work to ti productions followed 
by t — > £ + X (£ = e or /i) at hadron colliders, and investigate relevant non-standard 
top-quark couplings tig, tigg and tbW based on the effective Lagrangian which is 



constructed with SM particles alone. So far, a number of authors have studied 
this issue in similar frameworks to extract BSM signals at hadron colliders [5]-|25]. 
Since many of them are however focusing mainly on CP-violating observables, we 
study top-quark couplings which induce CP-conserving as well as CP-violating in- 
teractions. We first give updated constraints on the top-gluon couplings coming 
from the latest data on the tt total cross sections at the Tevatron and the LHC. 
Then we perform an optimal-observable analysis (OOA) to estimate the expected 
statistical uncertainties in measuring the top-gluon and top-W 7 couplings that con- 
tribute to pp — > tiX — > £ + X' under a linear approximation for current and future 
LHC experiments, which is our final goal here. 

2. Framework 

The effective-Lagrangian approach, in which its low-energy form reproduces the 
SM, is one of the most promising methods to describe new-physics phenomena 
when the energy of our experimental facility is not high enough to produce new 
particles. Assuming any non-SM particles too heavy to appear as real ones, the 
effective Lagrangian is given as 

i 

where Csm is the SM Lagrangian, Oi mean SU(3) x SU{2) x U{1) gauge-invariant 
operators of mass-dimension 6 involving only the SM fields fZE[ 127] and their coef- 
ficients Ci parameterize virtual effects of new particles at an energy less than the 
assumed new-physics scale /lo In this framework, all the form factors related to C,; 
are dealt with as constant parameters, without supposing any specific new-physics 
models. 

All the dimension-6 operators were initially listed by Buchmuller and Wyler |26j . 
However, it was pointed out that some operators there are related with others 
through equations of motion, which means that they are not independent of each 
other [28]. Then, the whole related operators were rearranged to get rid of this 
redundancy in Refs. (29J, [30] . 



f' 1 Since dimension-5 operators violate the lepton- number conservation, they are not treated 
hereafter. Therefore, we deal with dimension-6 operators, which give the largest contributions in 
relevant processes. 



Following the notation of Ref. [29], the relevant effective Lagrangian which 
describes non-standard interactions of the third generation for the parton-level 
process qq/gg — > tt — > £ + X is given in [23] as 
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where g s and g are the SU(3) and SU(2) coupling constants, Pl/r 
d v , d A and are form factors defined as 



(4) 

^ T 7 5 )/2, 



dy = 


gsA2 Ke{C uG(f> ), 


d A = 


g s A 2 lm ^«G^' 


fl = 


V tb + jp, 


fi = 


v 2 

2A 2 ' 




I, 2 

./of 33 * V 




./^ r 33 V 



(5) 



with v being the Higgs vacuum expectation value and Vtb being (tb) element of 
Kobayashi-Maskawa matrix. In particular, dy and d A are the so-called chromo- 
magnetic- and chromoelectric-dipole moments, respectively. 

In the following calculations, we use the above effective Lagrangian for top- 
quark interactions and the usual SM Lagrangian for the other interactions which 
are not affected by top quarks hereafter. 

3. Current bounds on dy and d A 

Experimental constraints on dy and d A were derived in our preceding works j20j [24] . 
Since LHC data of top-pair productions have been updated, however, we also re- 
derive new limits on dy and d A combining the recent data from the Tevatron for 
\fs = 1.96 TeV and the LHC for y/s = 7 TeV assuming m t = 172.5 GeV: 



a cxp = 7.50 ± 0.48 [pb], (CDF [31]) 



= 7.56±g : g [pb], (DO [32]) 

= 18711? [Pb], (ATLAS [33]) 

= 165.8 ± 13.3 [pb]. (CMS [M]) 

The corresponding theoretical predictions in the SM taking into account QCD 
higher-order corrections have been given by several authors [35] [35]. Instead of 
averaging them, we take here the following values: 

a SM = 7.291°,;^ [pb] for Tevatron, (6) 
= 1651$ [pb] for LHC, (7) 

which were used in [31], [33] . They are consistent with the others, and therefore 
this choice does not cause any serious problem to our results shown below except 
for some tiny difference. 

Adding those theoretical uncertainties to the above experimental errors quadrat- 
ically, we compare thus-obtained "effective" experimental data 

< = 7.50l°;| [Pb] for CDF, 

= 7 - 56 H:o2 [p°] for D °> 

= 187l| [pb] for ATLAS, 
= 1661^ [pb] for CMS, 

as our inputs with 

a = cr SM + Aa(d v ,d A ), (8) 

where ctsm represents the central values in the above equations (0[7]), Acr(dy, d A ) 
is the non-SM contribution calculated based on the effective Lagrangian, and the 
explicit forms of the corresponding parton-level quantities, i.e., Aa g g(d v , d A ) and 
Aa gg (dy, d A ) are found in the first article of Ref. [20]. 

Numerical analyses are performed the same manner (including the parton- 
distribution functions) as in [20] and the whole result is shown in Fig{TJ The 
shaded part there expresses the region of dy and d A allowed by the above Teva- 
tron and LHC data altogether. This shows clearly that both dy and d A are well 
restricted though there remains some small room for non-SM contribution. 
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Figure 1: The dy,A region allowed by the Tevatron and LHC data altogether (the 
shaded part). The solid curves, the dashed curves, the dash-dotted curves and the 
dash-dot-dotted curves are respectively from CDF, DO, ATLAS and CMS data. 



4. Optimal-observable analysis 

The optimal-observable analysis (OOA) is a method that could systematically es- 
timate the expected statistical uncertainties of measurable parameters. Here we 
do not go into its detailed description and just show how to know the uncertain- 
ties thereby, as more theoretical aspects of this method are presented in other 
papers [4"0]-[4"3]. 

First, in order to apply the OOA to the process pp — > tiX — > £ + X', we express 
its differential cross section (the angular and energy distribution of the charged- 
lepton £ + ) which we derived in our previous paper [23] as follows: 

— — — = fsu(Ei, cos 6^) + d v f dv (E e , cos 6^) + d R f dR (E e , cos 6^) 

dEpd cos &£ 

+ d 2 v f d 2 v (E e , cos9 e ) + d 2 A f d 2 A (E e , cos0<) + • • • , (9) 

where f SM (Ei, cos 9i) denotes the SM contribution, all the other fi(E(, cos 9i) de- 
scribe the non-SM terms corresponding to their coefficients, and d R is defined as 

d R = Re(f?)M w /m t . (10) 

The explicit forms of fi(E#, cos 9%) at the parton level are easily found in the relevant 
formulas in 



Since the magnitude of dy and d>A has been shown small as in FigJH we neglect 
any contribution from terms quadratic (or higher) in those non-SM parameters 
hereaftero Note that all d,A contributions disappear under this linear approxima- 
tion, because there is no term proportional to cIa- Thus, the angular and energy 
distributions of a decay lepton are written as 

dap 
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0i (cos 0*) + d v g 2 (cos 9p), (11) 
h(Ep) + d v h 2 {Ep) + d R h 3 (E e ), (12) 
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where gi(cos9p) and hi(Ep) are given by 

«(co.fc) = JdE t /,(«,,«.*), hi (E t ) =/dcosft ME t , mS e,) 

with i = 1,2 and 3 corresponding to I = SM, dy and d R , respectively. Here should 
be one comment about the angular distribution: We can thereby probe exclusively 
the dy term, since any contribution from da disappears within our approximation 
as a result of the decoupling theorem found in [15]-|4T]. 
Secondly, we calculate the following matrices: 
use I a a 9i(cos6p)g j (cos6p) 

My = / cicos^ — J — («, j = 1, 2), (13) 



g^cosOp) 



M « s /^'^fUr 1 (! ' j=1 ' 2 - 3) (14) 

and their inverse matrices , all of which are apparently symmetric. Then the 
statistical uncertainties for the measurements of couplings dy and d R could be 
estimated by 

\Sd v \ = yJXyjtjNt = y r Xjjl (15) 
through the angular distribution, and 



\5d v \ = \ Xgap/Np = JXi/L , (16) 



\5d R \ = ^Xgap/Np = yJXg/L (17) 

via the energy distribution, where ap, Np and L denote the total cross section, the 
number of events and the integrated luminosity for the process, respectively. 



" 2 The quadratic terms of fi-p have already been neglected in our previous paper [24]. Indeed, 
studies at the Tevatron suggest that those contributions are tiny: See [H] for the latest data. 



We are now ready to carry out necessary numerical calculations. Below we show 
the elements of M C ' E computed for y/s = 7, 8, 10, and 14 TeV, assuming m t = 173 
GevE 



(1) The angular distribution 

(1-1) ^ = 7 TeV 

M{ x = 23.102, M{ 2 = -245.412, M 2 C 2 = 2607.340. 
(1-2) v/i = 8 TeV 

M{ x = 33.234, M{ 2 = -353.598, M 2 C 2 = 3762.753. 
(1-3) v/i = 10 TeV 

M{ x = 59.333, M{ 2 = -632.179, M 2 C 2 = 6736.735. 
(1-4) = 14 TeV 

M c x \ = 134.052, M X2 = -1428.300, M 2 C 2 = 15220.286. 

(2) The energy distribution 

(2-1) v/i = 7TeV 

M X1 = 23.102, Mf 2 = -245.412, = 0.000, 

Mg = 2607.658, M| = -1.974, Mg = 15.323. 

(2-2) v/i = 8TeV 



M£ = 33.234, Mf 2 = -353.598, Mg = 0.000, 
Mg = 3763.252, M| = -2.880, Mg = 21.124. 



(2-3) v/i = 10 TeV 



Mg x = 59.333, Mg — —632.179, = 0.000, 

Mg = 6737.696, M| = -5.120, Mg = 35.233. 



" 3 We performed our analysis of the ti total cross section for m f = 172.5 in §3 since this value 
was assumed in the Tevatron and LHC data. Hereafter, however, we would like to use m f = 173 
GeV considering the latest world average |48j . 



(2-4) v/i = 14 TeV 

Mf x = 134.052, M E = -1428.300, Mf 3 = 0.000, 
M E = 15222.443, M E = -10.953, M E = 72.264. 

Here, all Mf 3 became zero for the same reason as vanishing da terms in the angu- 
lar distribution, i.e., the decoupling theorem [45]- (47]. Using the inverse matrices 
calculated from the above elements and Eq.(ll5|) - Eq. (ll7[) . we can estimate the 
statistical uncertainties of the relevant couplings. 

Before giving the results, however, we should comment about a stability problem 
in inverse-matrix computations: We noticed that the results fluctuate to a certain 
extent (beyond our expectation) depending on to which decimal places of M C,E we 
take into account as our input data. Therefore we calculate X C,E for the above 
M C,E and also for what are obtained by rounding those M C,E off to two decimal 
places. We then use their mean values to get 5dy and 5dn with "errors", which 
are the differences between the mean values and the maximum/minimum. Let us 
explain these treatments more specifically by taking (1-1) as an example: For those 
different inputs, we get 

^/Xf 2 = 1.73 for Ml x = 23.102, 
= 2.57 for M{ x = 23.10, 

which lead to the mean value 2.15, the maximum 2.57, and the minimum 1.73, and 
resultant errors ±0.42, which are derived by 2.57 — 2.15 and 1.73 — 2.15. These 
uncertainties are of course independent of the statistical ones which the OOA gives, 
i.e., Sdv,R themselves, so we should show both uncertainties separately as different 
errors. 

This instability is due to the large contribution from the dy coupling, which 
is dependent of gluon momenta and could therefore dominate. That is, we are 
consequently forced to treat large and small figures simultaneously in one matrix. 
In fact, we confirmed that any similar problem does not appear if we carry out 
computations for M E and X E without dy terms, in which there is no big difference 
between the SM- and dn-term contributions. 



We also have to explain how we can take into account QCD higher-order cor- 
rections: All the numerical computations of in this section were done with 
the tree-level formulas. In order to include QCD corrections there, we multiply 
the tree cross sections by the .fT-factor. This factor disappears in the combination 
X^- E ai and remains only in Nf, (= Lai) when we estimate 5dy y R- Therefore the 
luminosity L discussed in the following should be understood as an effective one 
including K (and also the lepton detection efficiency eg). 

We now show the whole results in Table [1] and Table |2] for the angular and 
energy distributions, respectively. The uncertainties we encountered here are hard 
to eliminate, but still the results will tell us the necessary luminosity for reaching the 
precision which we aim to realize. For example, we need at least L ~ 1500 pb _1 (= 
1.5 fb _1 ) in order to achieve \8dy\ ~ O(10 -2 ) in case of measuring the angular 
distribution at the LHC whose colliding energy is 8 TeV. As mentioned, those L 
should be divided by K ~ 1.5 and eg. If we assume q = 0.5, the resultant L 
increases slightly, which however hardly affects our conclusion. 



[TeV] 


\5dv \ x vT 


7 


2.15 ±0.42 


8 


2.25 ±0.96 


10 


1.12 ±0.12 


14 


0.73 ±0.02 



Table 1: Estimated statistical 
uncertainties of d v (multiplied by 
\/L) from the angular distribu- 
tion of a decay lepton. 



[TeV] 


\5dy\ x \/L 


\5dn\ x vT 


7 


1.86 ±0.28 


0.35 ±0.02 


8 


1.70 ±0.51 


0.32 ±0.05 


10 


0.98 ±0.08 


0.22 ±0.01 


14 


0.65 ±0.02 


0.15 



Table 2: Estimated statistical uncertainties of 
dy and d^ (multiplied by \/L) from the energy 
distribution of a decay lepton. 



Here we performed all the angular and energy integrations to compute M C,E 
over the full kinematically-allowed ranges for simplicity, but it will be useful to 
see how our results are affected by specific experimental conditions. As a typical 
example, we take the transverse-momentum cut for £ + in the energy distribution: 
Pn > p™™, which leads to 

-^l-{p™/E e y < cos6 e < +^/l-{p™ n /E i y 



with E t > pfi n . The results for p™ n =20 GeV (see, e.g., [331 M are 

\5d v \ (xVZ) = 2.30 ±0.51 and \5d R \ (xy/Z) = 0.69 ± 0.12 (18) 

for a/s =8 TeV, where our results with no experimental cuts were most unstable. 
We find there is no change in the instability of 8dy while that of 5dn has increased. 
However, the relative size of it compared to 6d R itself (= 0.69) is almost the same 
as that in the Table H 16 % -»■ 17 % (30 % -»■ 22 % in <Jd v ). Therefore, they tell 
us that introducing a lepton-pp cut will not seriously affect the stability of X C ' E , 
although the expected uncertainties themselves become a bit larger, which we can 
take into account through the detection efficiency as will be discussed in the last 
section. 

5. Summary 

Possible non-standard top-gluon and top-W couplings were studied for hadron- 
collider experiments in a model-independent way. Those couplings are derived as 
parameters which characterize the effects of dimension-6 effective operators based 
on a scenario of Buchmiiller and Wyler [26]. More specifically, we studied dy,A 
couplings which contribute to top-pair productions and d R coupling which affects 
leptonic decays of a top quark as non-standard interactions. 

First, we updated the bounds of dy and d& from the current Tevatron and 
LHC data of top-pair productions and showed that the allowed area of dyA be- 
came narrower around (dy,dA) = {0,0). Based on this observation, we adopted a 
linear approximation for the non-standard parameters in the subsequent optimal- 
observable analysis (OOA) of pp — > tiX — > £ + X'. In this approximation, all the 
contributions from rf^ were neglected because the leading order of d^ terms is 
quadratic. Therefore, the OOA was performed for dy and d R to estimate their 
expected statistical uncertainties using the angular and energy distributions of the 
final lepton. 

Measuring the energy distribution could study both dy and d R at the same time 
with a higher precision than the case of the angular distribution. Meanwhile, since 
the angular distribution is affected by only dy contributions within our approxi- 
mation, this one is suitable for exploring the mechanism of top-pair productions 



exclusively though the precision of dy is slightly lower than that from the energy 
distribution. 

We found that there were some ambiguities in inverse-matrix computations 
depending on where to round the input data (i.e., M^ E -elements) off, and con- 
sequently estimated statistical errors were affected thereby to a certain extent. 
However, our conclusion about the necessary luminosity and new-physics scale A 
hardly receive that serious influence, because those ambiguities do not change the 
order of the results. 

As for future prospects, if the LHC is going to be steadily upgraded and L = 
500 ftT 1 is achieved at y/s = 14 TeV, both dy and djt could be determined with 
\Sdy\ ~ O(10~ 3 ) and |<Mr| ~ O(10~ 4 ) from the energy distribution: Here, \Sdy\ ~ 
0.001 ~ 0.0001) means that contribution from effective operator Of^^ (C^fv) 

is suppressed by A > 7 (20) TeV if the center value of measured dy (da) is close 
to zero and C^q± (C^) ~ 1, which we estimated from Eqs.flH]) and ( TTUj) . 

Finally, we did not discuss here any realistic experimental conditions like effects 
of event-selection cuts and/or backgrounds except for the discussion in the last 
paragraph of the preceding section. This was partly because it was not our purpose 
to go into such details in this letter, although it is important since they decrease 
actual event number, but also because the OOA can work under any experimental 
conditions. All we have to do is to take account of them properly and accept the 
estimated luminosity as an effective one including all those effects. This is closely 
related to the detection efficiency ei, which we discussed at the end of the previous 
section. Indeed we are ready to re-calculate every quantity obtained in this letter, 
correspondingly to any lepton angle- and/or energy-cut conditions since we gave 
all the formulas used here analytically in Ref . [24J . 
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